Ground beetle (Carabidae) abundance data from north-eastern Australian rainforest, between June 2008 - January 2010, using pitfall traps

Understanding how the environment influences patterns of diversity is vital for effective conservation management, especially in a changing global climate. While assemblage structure and species richness patterns are often correlated with current environmental factors, historical influences may also be considerable, especially for taxa with poor dispersal abilities. Mountain-top regions throughout tropical rainforests can act as important refugia for taxa characterised by low dispersal capacities such as flightless ground beetles (Carabidae), an ecologically significant predatory group. We surveyed flightless ground beetles along elevational gradients in five different subregions within the Australian Wet Tropics World Heritage Area to investigate (1) whether the diversity and composition of flightless ground beetles are elevationally stratified, and, if so, (2) what environmental factors (other than elevation per se) are associated with these patterns. Generalised linear models and model averaging techniques were used to relate patterns of diversity to environmental factors. Unlike most taxonomic groups, flightless ground beetles increased in species richness and abundance with elevation. Additionally, each subregion consisted of distinct assemblages containing a high level of regional endemic species. Species richness was most strongly positively associated with the historical climatic conditions and negatively associated with severity of recent disturbance (treefalls) and current climatic conditions. Assemblage composition was associated with latitude and current and historical climatic conditions. Our results suggest that distributional patterns of flightless ground beetles are not only likely to be associated with factors that change with elevation (current climatic conditions), but also factors that are independent of elevation (recent disturbance and historical climatic conditions). Variation in historical vegetation stability explained both species richness and assemblage composition patterns, probably reflecting the significance of upland refugia at a geographic time scale. These findings are important for conservation management as upland habitats are under threat from climate change.

Secchi depth maps for the Baltic Sea, 2003-2012, with links to GeoTIFs

Secchi depth is a measure of water transparency. In the Baltic Sea region, Secchi depth maps are used to assess eutrophication and as input for habitat models. Due to their spatial and temporal coverage, satellite data would be the most suitable data source for such maps. But the Baltic Sea's optical properties are so different from the open ocean that globally calibrated standard models suffer from large errors. Regional predictive models that take the Baltic Sea's special optical properties into account are thus needed. This paper tests how accurately generalized linear models (GLMs) and generalized additive models (GAMs) with MODIS/Aqua and auxiliary data as inputs can predict Secchi depth at a regional scale. It uses cross-validation to test the prediction accuracy of hundreds of GAMs and GLMs with up to 5 input variables. A GAM with 3 input variables (chlorophyll a, remote sensing reflectance at 678 nm, and long-term mean salinity) made the most accurate predictions. Tested against field observations not used for model selection and calibration, the best model's mean absolute error (MAE) for daily predictions was 1.07 m (22%), more than 50% lower than for other publicly available Baltic Sea Secchi depth maps. The MAE for predicting monthly averages was 0.86 m (15%). Thus, the proposed model selection process was able to find a regional model with good prediction accuracy. It could be useful to find predictive models for environmental variables other than Secchi depth, using data from other satellite sensors, and for other regions where non-standard remote sensing models are needed for prediction and mapping. Annual and monthly mean Secchi depth maps for 2003-2012 come with this paper as Supplementary materials.

Responses of three tropical seagrass species to CO2 enrichment

Increased atmospheric carbon dioxide leads to ocean acidification and carbon dioxide (CO2) enrichment of seawater. Given the important ecological functions of seagrass meadows, understanding their responses to CO2 will be critical for the management of coastal ecosystems. This study examined the physiological responses of three tropical seagrasses to a range of seawater pCO2 levels in a laboratory. Cymodocea serrulata, Halodule uninervis and Thalassia hemprichii were exposed to four different pCO2 treatments (442-1204 µatm) for 2 weeks, approximating the range of end-of-century emission scenarios. Photosynthetic responses were quantified using optode-based oxygen flux measurements. Across all three species, net productivity and energetic surplus (PG:R) significantly increased with a rise in pCO2 (linear models, P < 0.05). Photosynthesis-irradiance curve-derived photosynthetic parameters-maximum photosynthetic rates (P max) and efficiency (alpha) also increased as pCO2 increased (linear models, P < 0.05). The response for productivity measures was similar across species, i.e. similar slopes in linear models. A decrease in compensation light requirement (Ec) with increasing pCO2 was evident in C. serrulata and H. uninervis, but not in T. hemprichii. Despite higher productivity with pCO2 enrichment, leaf growth rates in C. serrulata did not increase, while those in H. uninervis and T. hemprichii significantly increased with increasing pCO2 levels. While seagrasses can be carbon-limited and productivity can respond positively to CO2 enrichment, varying carbon allocation strategies amongst species suggest differential growth response between species. Thus, future increase in seawater CO2 concentration may lead to an overall increase in seagrass biomass and productivity, as well as community changes in seagrass meadows.

Modelled spatial and seasonal distribution of Blue Whiting (Micromesistius poutassou) larvae in the North-East Atlantic (1951 to 2005)

Blue whiting (Micromesistius poutassou, http://www.marinespecies.org/aphia.php?p=taxdetails&id=126439) is a small mesopelagic planktivorous gadoid found throughout the North-East Atlantic. This data contains the results of a model-based analysis of larvae captured by the Continuous Plankton Recorder (CPR) during the period 1951-2005. The observations are analysed using Generalised Additive Models (GAMs) of the the spatial, seasonal and interannual variation in the occurrence of larvae. The best fitting model is chosen using the Aikaike Information Criteria (AIC). The probability of occurrence in the continous plankton recorder is then normalised and converted to a probability distribution function in space (UTM projection Zone 28) and season (day of year). The best fitting model splits the distribution into two separate spawning grounds north and south of a dividing line at 53 N. The probability distribution is therefore normalised in these two regions (ie the space-time integral over each of the two regions is 1). The modelled outputs are on a UTM Zone 28 grid: however, for convenience, the latitude ("lat") and longitude ("lon") of each of these grid points are also included as a variable in the NetCDF file. The assignment of each grid point to either the Northern or Southern component (defined here as north/south of 53 N), is also included as a further variable ("component"). Finally, the day of year ("doy") is stored as the number of days elapsed from and included January 1 (ie doy=1 on January 1) - the year is thereafter divided into 180 grid points.

n-alkane carbon isotope ratios in plants and recent lake sediments from Cameroon

Trees and shrubs in tropical Africa use the C3 cycle as a carbon fixation pathway during photosynthesis, while grasses and sedges mostly use the C4 cycle. Leaf-wax lipids from sedimentary archives such as the long-chain n-alkanes (e.g., n-C27 to n-C33) inherit carbon isotope ratios that are representative of the carbon fixation pathway. Therefore, n-alkane d13C values are often used to reconstruct past C3/C4 composition of vegetation, assuming that the relative proportions of C3 and C4 leaf waxes reflect the relative proportions of C3 and C4 plants. We have compared the d13C values of n-alkanes from modern C3 and C4 plants with previously published values from recent lake sediments and provide a framework for estimating the fractional contribution (areal-based) of C3 vegetation cover (fC3) represented by these sedimentary archives. Samples were collected in Cameroon, across a latitudinal transect that accommodates a wide range of climate zones and vegetation types, as reflected in the progressive northward replacement of C3-dominated rain forest by C4-dominated savanna. The C3 plants analysed were characterised by substantially higher abundances of n-C29 alkanes and by substantially lower abundances of n-C33 alkanes than the C4 plants. Furthermore, the sedimentary d13C values of n-C29 and n-C31 alkanes from recent lake sediments in Cameroon (-37.4 per mil to -26.5 per mil) were generally within the range of d13C values for C3 plants, even when from sites where C4 plants dominated the catchment vegetation. In such cases simple linear mixing models fail to accurately reconstruct the relative proportions of C3 and C4 vegetation cover when using the d13C values of sedimentary n-alkanes, overestimating the proportion of C3 vegetation, likely as a consequence of the differences in plant wax production, preservation, transport, and/or deposition between C3 and C4 plants. We therefore tested a set of non-linear binary mixing models using d13C values from both C3 and C4 vegetation as end-members. The non-linear models included a sigmoid function (sine-squared) that describes small variations in the fC3 values as the minimum and maximum d13C values are approached, and a hyperbolic function that takes into account the differences between C3 and C4 plants discussed above. Model fitting and the estimation of uncertainties were completed using the Monte Carlo algorithm and can be improved by future data addition. Models that provided the best fit with the observed d13C values of sedimentary n-alkanes were either hyperbolic functions or a combination of hyperbolic and sine-squared functions. Such non-linear models may be used to convert d13C measurements on sedimentary n-alkanes directly into reconstructions of C3 vegetation cover.